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ERROR ANALYSES OF RESONANT 
ORBITS FOR GEODESY 


B. C, Douglas 

C. F. Martin 
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R.G. Williamson 


ABSTRACT 

An error analysis of resonant orbits has indicated that attempts 
to recover high order resonant geopotential coefficients willbe seri- 
ously hampered by errors in the non-resonant geopotential terms. 
This effect, plus the very high correlations (up to .999) of the reso- 
nant coefficients with each other and the orbital period in single 
satellite resonant solutions, makes individual resonant orbits nearly 
valuless for geodesy. Multiple -satellite, single-plane solutions are 
only a slight improvement over the single satellite case. Independent 
determinations of resonant terms from low altitude satellites require 
multiple orbital planes and small beat -periods to reduce correlations 
and effects of errors of non-resonant geopotential terms. Also, these 
unmodeled parameter effects on low-altitude resonant satellites make 
the use of tracking arcs exceeding two to three weeks of doubtful 
validity. Because high-altitude resonant orbits are less affected by 
non-resonant terms in the geopotential, longer tracking arcs can be 
used for them. 




This project was sponsored by NASA under 
Contract No. NAS 5-9756-135. 





ERROR ANALYSES OF RESONANT 
ORBITS FOR GEODESY 


1. INTRODUCTION 

This study concerns determining the realistic accuracy to which certain 
geopotential coefficients can be recovered by observing satellites on resonant 
orbits. The results will be useful to investigators contemplating either resonant 
satellite geodetic missions or analysis of existing resonant satellite orbits. 

Reference 1 presented a "noise only" error analysis of resonant orbits. 
Designed to yield only relative information, that study did not produce absolute 
predictions of the quality of determinations of gravity coefficients from resonant 
orbits. Error analyses are almost always very optimistic unless we include the 
effects of errors in the unmodeled or unadjusted parameters of the problem. 

The resonant satellite geodesy problem is an extreme example of this phenome- 
non because the "noise only" standard deviations are as much as five orders of 
magnitude optimistic in some cases. 

This study includes effects of errors in the following parameters : 

• non^resonant geopotential coefficients; 

• station location; 

• unadjusted resonant geopotential coefficients; and 

• GM, the product of the Earth's mass and the gravitational constant. 

Errors in all of these quantities have large effects on the solution. Indeed, for 
low altitude (or high order) resonant satellites in a single plane, errors in the 
non-resonant geopotential coefficients entirely dominate the solution. 

Finally, the causes of the large effects of the errors in the unadjusted 
parameters have been sought in terms of the physical characteristics of resonant 
orbits. 

2. SHALLOW RESONANCE 

For several reasons this study is confined to orbits in "'shallow" resonance - 
approximate rather than exact or "deep" resonant orbits. First, very deep 



resonance is not required to achieve large perturbations such as several kilom- 
eters along track; second, deeply resonant orbits are achievable only after several 
orbit corrections. From a mission standpoint, such a requirement me mis greatly 
increased spacecraft complexity. Finally, this study lias shown that for low mid 
medium altitude satellites deeply resonant orbits are less desirable for satel- 
lite geodesy than orbits in shallow resonance because they are affected more by 
unmodeled parameters. 


The terms “shallow" and "deep" resonance can be given a reasonably pre- 
cise meaning in terms of satellite beat period. The condition for an exactly 
repeating ground track on the average, or deep resonance, is 


D 


t M 


w 



<0 


0 , 


where the subscript c denotes secular rate, and s is the integer number of 
revs/day made by the satellite. <D is the drift rate of the orbit, G and o are 
the rates of rotation of perigee, the node and the earth. M is the mean motion 
of the satellite.) (See Reference 1.) For shallow resonance, D is not zero but is 
small compared to the satellite mean motion. For example, the GEOS-II satel- 
lite orbiting the Earth approximately thirteen times per day has 


D rr „ 57. 2°. day, 


compared with 


M s 4619° day. 


The beat period, 360 °/D, is a measure of the time for the ground track to re- 
peat. For GEOS-II this is 6'.'3. 

The most important parameters to consider when selecting orbits for reso- 
nant satellite geodesy are beat period, inclination, and eccentricity. The in- 
fluence of these parameters will be demonstrated. 

The solutions for perturbations of orbital elements for deeply resonant 
orbits are different from those for non-resonant orbits. The non-resonant 
solutions assume small perturbations, and we can treat them as linear forced 
oscillations about an intermediary, or reference, trajectory. Of course, deep 
resonance violates tins condition because the satellite may be perturbed up to 
90° along track, not a small change. However, Palmiter and Gedcon in Reference 
2 have shown that the non-resonant perturbation solutions are quite accurate for 
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small beat periods and are in error only a lew percent lor beat periods near 
100 days. Thus* we define shallow resonance as corresponding to repeating 
ground track lx ‘at periods of about 100 days or less and reserve the term deep 
resonance lor longer beat periods. However, using non-resonant formulas on 
orbits with beat periods near 100 days requires considerable caution because 
the observable resonant perturbations on an orbit may vary in magnitude by a 
factor of 100 or more. Obviously, even a one percent error in calculating the 
larger effects may seriously degrade the quality of determining the sources 
of the smaller effects. Numerical integration is probably the only reasonable 
procedure for beat periods greater than approximately fifty days. 

The graphs in Reference 1 show how the acceleration of mean anomaly 
caused by dominant resonant geopotential coefficients varies with inclination. 
These graphs indicate that high inclinations, i . 40°, are required to achieve 
large perturbations for medium and low altitude resonant orbits, s> 8 revs/day. 
Our present report exhibits similar information for several selected representa- 
tive cases, yet shows variation with eccentricity and bent period as well. To 
accomplish this, we lumped the along -track effects of many or all of the resonant 
coefficients into a Root-Sum-or-Squarcs (RSS). The along-track effect was 
estimated from 

/(along* track) „ n(A> * AM t //boos i) 
where A . . AM , and AQ appear in Reference 3 as 
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The quantity is the integral of Sf. with respect to its argument. The 

p,q indices identify particular harmonic ’components of a spherical harmonic 
The functions Jy mp (i) and Gy (e) are defined in Reference 3. Gp pq (e) 
is 0( o I n i ) for low and moderate e . 

In Reference X and elsewhere, the resonant harmonic components (t, m , p , 
q ) are given by the condition that 


4 - 2p t q 


m 

ft 


For example, the resonant ( £, rn , p , q) sets for a satellite with mean motion at 
or near fifteen teims per day are 
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15, 15, (i, -2 
15, 15, 7, 0 

15, 15, B, 2 
# » • • 

» » i ♦ 

• * * « 

Hi, 13, 7, -1 
1G, 15, 8, 1 

i • I » 

• » t i 

• f I I 

In the present work, all resonant (i , m , p , q) sets to (i , m) - (30, 30) were 
computed and included in the perturbation ealouation^. 

As a consequence of the treatment of tea serai harmonic perturbations as 
linear forced oscillations we also speak of a beat period for each resonant (4, 
m , p , q) component calculated from 360°/J\ , where 

f -mp q 

^ mpq ^ “ ^P) a c 1 ^ "* ^P + C 1)M C 1 m(H c - V ) • 

For a circular orbit, or one at the critical inclination, all resonant terms have 
the same beat period, the repeating groundtrack value. This does not apply to 
eccentric orbits at arbitrary inclination bvie? use of the secular advance of 
perigee caused by J 2 . The resonant (4,m , sets for the example above 
would yield the corresponding drift rates D^ npq : 

^15 , 15 . 6 , ~2 ’ ‘^'e + K + ** (H c ~ & ) * 

^15, IS* 7,0 *«e+Mc+ 1S C^-&) 1 

■^ 16 , 15 , 7,-1 f W c 1 15(fi c - 60, etc, 
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Only in the q 0 circular case is the beat period of a resonant (?, m, p ,q) har- 
monic component the same as the repeating groundtrack beat period. The drift 
rates and the resulting beat periods for the various components vary according 
to the contribution of ( f ' - 2p ) ‘ . Of course, the larger the beat period, the 
greater the effect oi any particular ( i’, m, p,q) component. For a beat period 
> 5 days the term with the quadratic divisor in A Mf;, l)pq will dominate; thus, the 
effect will increase quadratically with beat period as long as Equations 2-1 are 
valid. 


is always small compared to M e or o ; thus, whether (•!- 2p) is 0, 2, or 
3 does not have much effect for very shallow resonance (< 10 days). However, 
for somewhat deeper resonance, especially for beat periods >30 days, is 
comparable to (4 - 2p + q)M c + m(U c - 0). The presence or absence of may 
drastically alter the beat period of a particular (i, m , p ,q) component compared 
4 to the repeating groundtrack beat period. Figure 2-1 show how the RMS of the 
beat periods of the resonant terms varies with repeating ground track beat 
period, and inclination for an orbit resonant with 10th order of geopotential terms, 
The identifying numbers on each curve indicate the repeating groundtrack beat 
period for the case. 



Figure 2- 1 — Effect of Inclination on RMS Beat Period of Resonant ("t, m, p, q) Components 

for S * REVS/DAY 
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In Figure 2-1, first note that the ten day repeating groundtraek beat period 
case shows almost no variation of resonant beat periods with inclination. Longer 
repeating groundtraek beat periods show the RMS beat period varying consider- 
ably both wit nelination and beat period except at i ™ (33?45, the critical in- 
clination. At i - (>3?45, ' 0, and all beat periods are the same. Of course, 

all resonant terms beating with the same frequency is not conducive to separat- 
ing the terms from one another. For high altitude satellites, this situation is ag- 
gravated at all inclinations because * c declines rapidly with an increasing semi- 
major axis. 

Finally, an orbit with a large negative beat period is close in semi-major 
axis to an orbit with a large positive beat period, Thus, the values of perturba- 
tions vary irregularly at long beat periods because a small change in repeating 
groundtraek beat period may cause a large number of resonant (f., m , p , q) 
components to radically alter their beat period, even changing sign. Several fol- 
lowing graphs will depict this irregularity. 

Tlie effect of eccentricity on the resonance phenomenon is considerable. 
Reference 1 conveys that maximizing eccentricity maximizes resonant perturba- 
tions; however, an excessive eccentricity might be poor because too many terms 
might have perceptible effects. Table 2-1 fx'om Reference 4 shows an extreme 
example of this phenomenon. This is an analysis of some of the resonant per- 
turbations on 19G6 92a, a Russian twelve-hour communication satellite. Second 
order terms from (2,2) to (19,2) produce effects in the transverse (along-track) 
direction larger than 10 km! 


Figure 2-2 shows how the RSS of the dominant resonant perturbations along- 
track vary as a function of beat period for several eccentricities and inclinations 
for an orbit resonant with 10th order terms. The y axis shows the Root-Sum- 
of-Squares of the along-track perturbations of (10,10) through (15,10). 

A 10th order resonance was the basis of Figure 2-2; however, the along- 
track displacements are typical of orbits in shallow l'esonance. For high inclina- 
tion orbits of moderate to high eccentricity, a beat period of twenty days 
virtually guarantees effects of several km along-track. The effects for negative 
beat periods - mean motion less than an integer times the rotation rate of the 
Earth - are comparable. Thus, deep resonance is not required to produce easily 
observable effects and is even detrimental to determining geopotential coeffi- 
cients, as will be seen later. Figures 2-3, 2-4 and 2-5 show similar information 
for satellites resonant with fourth order terms. Note the very large effects for 
e = .55. 
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Table 2-1 


SATELLITE 66 92A 

S * 2 REV. /DAY PERIGEE HEIGHT * 269. KM. 
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REPEATING GROUND TRACK BEAT PERIOD, (DAYS) 

Figure 2-2-RSS Along-Track Perturbations for 
Resonant Terms with S - 10 REVS/DAY 





Figure 2-3-Effect of Inclination on Along-Track 
RSS Perturbations for S = 4 REVS/DAY, e = ,55 



RSS AMPLITUDE (METERS) 



BEAT PERIOD (DAYS) 

Figure 2-4-Effect of Inclination on Along-Track 
RSS Perturbations for S ”4 REVS/DAY, e - .25 
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Choosing orbital eccentricity therefore involves two important factors. 

First, the eccentricity must be large enough so that the eccentric resonant terms, 
l even for s odd and vice versa, have observable effects. However, if eccen- 
tricity is too large, a solution for any subset of the resulting numerous resonant 
terms will be corrupted by the unmodeled ones. Figures 2-G and 2-7 show the 
ratios of the RSS along-track effects due to circular and eccentric resonant 
terms for the tenth order resonant terms with £ u 15 for e ~ .05 and .1. Note 
that for e =- .1, the circular and eccentric resonant terms are comparable in 
importance. Large eccentricity is not required to make the eccentric resonant 
terms visible. Figure 2-8 for the fourth order resonance shows a similar re- 
sult. Note in Figure 2-8 the lack of irregularity in the curves. For a satellite 
resonant with fourth order terms, is very small and has little effect on beat 
periods. 

Figures 2-9 and 2-10 show how eccentricity affects the number of resonant 
terms that have lai'ge effects, These figures are plots of the ratio of the RSS of 
the resonant terms with t ~ m l 5 to resonant terms with 5 £ - m < 20 for an 

orbit x’esonant with tenth order terms. Note that for e - .2 and medium inclina- 
tion, i 2: 40°, this ratio is about 1:1. This is a poor situation. The high degree 
terms would certainly corrupt the determination of the terms with t 2 5. How- 
ever, at high inclinations, the terms with small (4-m) dominate, and the ratio 
increases to about ten for e ~ ,2 and 20-30 for e - .1. Thus, an eccentricity of 
about .1 appears to be a good choice for low altitude resonant orbits. It is high 



Figure 2*6~Effect of Inclination on the Relative Strength of Eccentric Resonant Terms; 

S - 10 REVS/DAY, e s. ,05 
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Figure 2*7- Effect of Inclination on the Relative Strength of Eccentric Resonant Terms; 

S « 10 REVS 'DAY, e .1 



BEAT PERIOD, (DAYS) 


Figure 2-8-Effect of Eccentricity and Inclination on the Relative Strength of Eccentric 

Resonant Terms for S - 4 REVS T) AY 


enough so that eccentric orbit resonant terms have effects as large as the 
circular orbit resonant terms, yet not so large that an unmanageable number of 
terms exist. A similar situation also generally applies to all other eccentric 
resonant orbits. 









BEAT PERIOD (DAYS) 


Figure 2-10— Effect of Inclination on the number of Resonant Terms with Large 
Perturbations for S s 10 REVS/DAY, e sr ,1 







We are not yet tempted to draw conclusions about the ideal orbit elements 
for resonant satellite geodesy, The conclusions would in many eases be quite 
erroneous. We ean only draw conclusions about the "best” orbit« after consider- 
ing the effect of the other important factor in the problem* namely the inevitable 
errors in the parameters not part of the solution vector. 


0. ERROR ANALYSIS OF RESONANT ORBITS 

Resonant orbits are extremely attractive for geodesy. Resonance enhances 
perturbations caused by tesseral harmonics and is especially important for high 
order terms. In the non-resonant case, terms beyond the eighth order may have 
an effect on an orbit of only a few meters. Furthermore, it is very difficult at 
present to extract spherical harmonic coefficients of this order from gravimetry, 
Thus, accurate coefficients <rom resonant orbits will provide important standards 
for evaluation of gravimetric results in addition to improving our capability to 
determine orbits. 


We have attempted to discover the degree of accuracy that we may deter- 
mine geopotential coefficients by using high quality observations of satellites in 
resonant orbits. As in many other "real world" analyses, the data quality is 
not the limiting factor; the errors in the unmodeled (unadjusted) parameters 
cause severe degradation of the solution over the "noise only" situation, A de- 
tailed description of how the effects of errors in unmodeled parameters are 
propagated into the estimates of the variances of the adjusted parameters appears 
in the Appendix, 

In this study, the adjusted parameters were generally a subset of the infinite 
number of resonant geopotential coefficients for the orbit in question, and the 
elements of the orbit or orbits used in the solution. The unmodeled parameters 
were 

• unadjusted resonant terms, 

• unadjusted non-resonant terms, 

• station location errors, and 

• uncertainty in GM. 

Table 3-1 lists the a priori uncertainties in these quantities and the tracking- 
station complement and schedule. 

A priori estimates of the errors of the unmodeled parameters are somewhat 
controversial because we tend to question the raw variances that may result from 
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Table 3-1 

Unadjusted Parameters 


S - 10 revs/day 

S - 4 revs/day 

All Cases 

Parameter 

A Priori 

Parameter 

A Priori » 

Parameter 

A Priori ^ 

C 

W H, in 

2*10 lf ’ 

C..« 

4 <10 11 

A(GM) 

10 

S H. in 

2 jf 10"" 18 

S. l4 

2.2x10 10 

ACOM *X 

-20m 

OlH, 10 

8x 10 19 

^9,4 

4 X10 11 

ACOM Y 

-20m 

® 1 J . 10 

gxlO' 1 '' 

S P,4 

1,0x10“" 

ACOM Z 

-20m 

^10,1) 

6xl0- 19 

c 

10,4 

s 

10,4 

2 xlO’* 11 

Ageopotcntial 

.25 (API.-SAO) 

S 16, 10 

(jxlO- 19 

0.5x10“" 

(non- 
res onant) 



Tracking Stations 


Wallops Station 


Latitude 

37?859 


Longitude (E.) 

Wrmrammmir: 

284.490 


Height above Spheroid (m) 


30.40 


Winkfield (UK) 
Hawaii 
Tanancrive 
Western Test Range 
Carnarvon 


51.445 

21.521 


-19.009 


37,500 


-24.903 


359.302 

202.003 

47.300 

237-502 

113,716 


70.00 

368.31 

1355.87 

195.03 

10,54 


Data quality; o range « 1,4m, n angles ~ 20 H @ 1 obs/sec. 


Tracking Schedule** 


Start 

Stop 

0“ 

4 d 

6 11 

91? 5 

12 d 

10 d 

18 d 

22 d 

24 d 

28 d 


30 d 




a least squares fit. Our a priori estimates of uncertainty represent what we 
believe is reasonable or will be attained in the near futoe. 

The a priori values of the unadjusted tenth order resonant terms have been 
taken as the expected order of magnitude of these terms. The earth’s gravity 
coefficients generally follow the rule 





0 



I 


where the overbars indicate normalised coefficients. For the unadjusted fourth 
order terms, essentially the same rule was used. This may be pessimistic be- 
cause there are published values of many fourth order coefficients. In both 
cases the list of unadjusted resonant terms was truncated at a degree for which 
't - m > 7. For low and moderate eccentricities, terms for which i - m > 7 have 
enormously reduced effects compared to terms for which £ -m or l = m + 1. 

Determining the error in the non-resonant geopotential terms is difficult. 

We chop ) to assume that the term by term differences between gravity models 
are some measure of their uncertainty. We have used e (APL 3.5 - SAO Ml), 
where e is a scale factor. Because the SAO Ml model is more recent and yields 
better results for orbit determination than the APL 3.5, e was taken as .25. The 
uncertainty in GM was taken as 1.T0 6 . 

We chose six tracking stations having good distribution in longitude and highly 
accurate tracking equipment. We have assumed thau GSFC Laser or FPQ-6 radar 
quality ranges were taken with cr ~ 1*4 m at a data rate of 1 obs/sec. The 

range 

angular data were assumed good to 20”; however, in the two to four week arcs 
considered, angular data of this quality have an entirely negligible effect on the 
solution. The a priori estimates of uncertainty of the station positions were 
taken as -20m in each center-of -mass coordinate. Biases in the range data were 
not considered because biases for these instruments are known to be much 
smaller than the station location error figure of -20 m . 

A few simulations of the recovery of resonant geopotential coefficients 
showed that data quality is not a limiting factor for resonant satellite geodery, 
at L ast within reasonable limits. The contribution of the pure noise to the 
estimates after two weeks of Lacking is always two to three orders of magnitude 
less than the order of magnitude of the coefficients themselves. Thus we did not 
evaluate factors such as the effect of adding more tracking stations and different 
data types. Resonance effects are so large that with modern tracking systems, 
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the random noise does not significantly affect the estimates of the coefficients 
compared to the effects of the* unmodeled parameters. 

Determining geopotential coefficients is plagued by the problem of separa- 
bility. The principal non-resonant uffv ot of a geopotential term of order m is 
an m times daily oscillation. Thus, wo require many different orbital pih'vss to 
provide die information for separating the effects of the various terms of a given 
order. 


A similar situation exists for circular resonant orbits or for resonant orbits 
at the critical inclination, wherein all resonant terms produce perturbations of 
the same frequency. Using one of these orbits, we cannot determine more than 
a single "lumped" coefficient. Contrarily, the frequency content of an eccentric 
resonant orbit, not at i i , with a boat period greater than ten to fifteen days, 
is considerable. The present authors and others have suggested the possibility 
of recovering a large number of geopotential coefficients from a single eccentric 
resonant orbit. However, simulated satellite tracking on a .1 eccentricity orbit 
having a -30 day beat period resonant with tenth order terms for four weeks has 
yielded very poor results due to errors in the unmodeled parameters, especially 
errors in the non-resonant geopotential coefficients. No significant information 
concerning individual resonant terms was obtainable from single resonant orbits. 


Failure to obtain useful information from a single eccentric resonant orbit 
led to attempted multiple-satellite solutions. Figure 3-1 compares the single- 
satellite solution with a three -satellite solution. The solution vector contains 
the orbital elements and the resonant terms from (10, 10) - (13, 10). The solu- 
tions shown for (10, 10) and (13, 10) are typical. The contribution to the solution 
uncertainty of the noise is off the bottom of the scale. 


The units on the y-axis of Figure 3-1 and of subsequent similar figures are 
those of standard deviations (o). On the left of each In? indicated by arrows are 
the a priori rj , and the or of the solution including the effects of the unmodeled pa- 
rameters. On the right side of the bars are contribution levels of the various un- 
modeled parameters. The total effect of all of these is their Root-Sum-of-Squares. 
Rather than give the effect of the error of each center -of-mass (COM) component, 
the contributions in each direction are combined as COM X, COM Y, COM Z. 

In Figure 3-1 we see the startling result that the o’ s of the fits are in some 
cases larger than the a priori values. But the o-’s shown include the effects of the 
unadjusted parameters. The noise -only o-’s are always two to three orders of 
magnitude smaller than the a priori information. The appearance of o’ s that 
include unadjusted parameter effects larger than the a priori values implies 
that the corrections from an attempted fit would be so bad that one wcii/ci have 
been better off not making the fit. For an example, if one knew a station position 
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within 20 meters and tried to improve the location by using near-Earth satellite 
data yet ignoring the effect of J 2 , he would obviously degrade his knowledge with 
such a fit. 

Figure 3-1 shows that although the noise-only solutions may be a great 
improvement over the a priori variances of the constants, including the effects of 
unadjusted parameters would not permit much more than order-of -magnitude in- 
formation for more than a few coefficients even with three satellites. Further- 
more, the correlations of the constants with each other and/or the orbital semi- 
major axes are in many cases .95 or higher. 

The unmodeled parameters causing the largest degradation of the solution 
are the unadjusted non-resonant terms. Their contribution to the total is de- 
noted ’’APL-SAO gravity.” 


The degree of the effect of errors in the non-resonant terms was unexpected; 
however, further investigations revealed the reason for this large effect. 

We know that errors in the gravity field prevent accurately determining 
orbital energy, thus along -track position. The resonance effect manifests itself 
almost entirely along-traek and is extraordinarily sensitive to semi-major axis 
for beat periods >10-15 days (as discussed in Section 2). This fact is also re- 
vealed by very high correlations > .9 between many of the » ^sonant constants and 
semi -major axes. Thus, we expect the errors of the non-resonant geopotential 
terms to have a large and increasing effect with increasing beat period. Figure 
3-2 illustrates this. Here we are comparing identical two-satellite solutions for 
tenth order resonant terms. In both cases, orbital eccentricities were 1/10, 
and inclinations were 90° during four weeks of simulated tracking. The long beat 
period multiple satellite solution with one beat period > 100 days has superior 
noise -only results to the short beat period case of ±10 days. However, the ef- 
fect of the unadjusted non-resonant gravity terms is vastly decreased in the 
dominant (10, 10) resonant term for the short beat period, seen by comparing 
the values denoted by ’’APL-SAO gravity.” Furthermore, as expected, the 
correlation coefficients for the resonant constants and the semi-major axes of 
the orbits io reduced for the small periods, some of them approaching .9 for 
the ±10 day beat period in contrast to about .97 for the long beat period case. 

Figures 3-1 and 3-2 show the real problem confronting the resonant satel- 
lite geodicist is finding orbital situations for combined solutions that minimize 
the effect of the unmodeled parameters. Figures 3-3, 3-4, and 3-5 show the 
results of attempts at reduction of unmodeled parameter effects by variation of 
various orbit and other parameters. 
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Figure 3-3 compares the results of two, two-satellite solutions. Eccentric- 
ities for case 1 were the same at .15 and for case 2 were set at .1 and .15. Two 
weeks of tracking were assumed for these polar orbits. A combined solution 
with different eccentricities sometimes produces improved results although 
there is no meaningful information with this change for terms with (E -m ) ;> 1, 

Another attempt with varied orbital geometry used identical eccentricity of 
.15 and inclination of 90°; however, arguments of perigee and node were 90° and 
135° apart, respectively. Again, varying the orbital elements produced a slight 
improvement only for small ( i - m ). 

Figure 3-4 shows the results of observing multiple satellites in the same 
plane for two and four week periods. Note that the value of a(C 10i 10 ) is worse 
for the longer arc. This phenomenon is often observed in orbit determination 
problems. Two long a tracking arc may degrade a solution because of the effects 
of the unmodeled parameters. 

Inclination significantly affects the resonance phenomenon. Figure 3-5 
shows the results of two-satellite solutions at 45° and 85°. The solution at the 
high inclination is superior for the dominant C 10 10 term, as would be expected. 
At the lower inclination, the effects of the various resonant terms are more 
nearly equal in magnitude than at the high inclination where terms with ( t -m) 
equal to 0 or 1 dominate. At the low inclination the effects of the unmode led 
resonant terms are very large for this reason. 

The single plane cases in Figures 3-1 through 3-5 have a very common 
factor: The errors in the non-resonant gravity terms dominate the estimates. 

We can control the effect of the unmodeled resonant terms by selecting eccen- 
tricity. Station location errors and GM error are seldom serious. The large 
effect of the unmodeled non-resonant terms is caused by the sensitivity of the 
resonance phenomenon to the orbital energy or period. 

Having failed to achieve meaningful results with a single plane, multiple 
satellite solution, we next attempted multiple orbital planes. Figure 3-6 com- 
pares the previous three -satellite, single-plane solution with a two -satellite, 
two-plane solution. The effects of the unadjusted non-resonant gravity terms 
have decreased greatly over the single plane solution, particularly for large 
(4 -m )- This multiple plane solution has been the first to yield substantial in- 
formation about a number of resonant coefficients. Even more important, the 
correlations of the geopotential coefficients with the orbital semi-major axes 
were reduced to .8 at the most. 

In the multiple plane solution we did not expect but did obtain an increased 
effect of the unadjusted resonant terms. However, one of the planes was inclinded 


26 







45®. Recall that low and middle inclinations are less favorable in this respect 
because the resonant perturbations in these orbits are more nearly equal for 
both small and large values of ( f - m ). 

The problem of correlations appears to be worse for high altitude orbits. 

A correlation coefficient of .999 was observed between C A 4 and C 7 4 for case 
one of Figure 3-7. furthermore, the correlation coefficients with the semi- 
major axis, were ,9 in some cases, in spite of the multiple orbit planes. There- 
fore, constants determined from this situation would have no meaningful appli- 
cability other than to these orbits. Many orbit planes are necessary for high 
altitude i*esonant orbits to yield independent estimates of the gravity field. 


1, CONCLUSIONS 

We have ascertained that the random noise level required for resonant 
satellite geodesy is not stringent; moreover, resonant satellite geodesy docs 
not require continuous tracking by numerous stations. However, data biases, 
like station locations may be important if excessive. Further, tracking arcs 
longer than a few weeks for low altitude resonant orbits are undesirable be- 
cause of the effects of errors in the non-resonant geopotential coefficients. 

We have also demonstrated that single, eccentric resonant orbits alone are 
not valuable to geodesy if more than order -of- magnitude estimates of geopoten- 
tial coefficients are required. The estimates of the errors of the geopotential 
coefficients obtained from a single resonant satellite are almost perfectly cor- 
related with each other as well as with the estimates of the error of the semi- 
major axis. Thus, the coefficients from a single satellite will be applicable only 
to the orbit from which they were obtained. This fact, while unfortunate for 
geodesy, suggests that for operational purposes, when resonances exist, we can 
obtain good results for prediction purposes by solving for one or two "lumped" 
resonant coefficients. We have in addition shown that the single plane multiple 
satellite solutions are only a slight improvement over single satellite solutions. 
As in the single satellite case, the controlling factors are the errors in the un- 
modeled non-resonant terms and the very high correlations of the coefficients 
with each other and with orbital semi -major axes. 

Finally, small (2:10 day) beat periods are desirable to uncouple the geopo- 
tential coefficients from the orbi'al semi-major axes. The low beat periods 
also mean that period errors of at least several minutes can be tolerated in a 
resonant satellite geodesy mission. 
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CASE 1 CASE 2 


Figure 3*7~Two Plane Solutions for 4^ Order Terms 
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Thus, resonant satellites in a single plane must be used in a combined 
solution with non-resonant satellites or gravimetry or both to make a meaning- 
ful contribution to geodesy Satellites of low to moderate eccentricity in several 
orbital planes (with some of high inclination), all resonant with terms of the 
same order, can produce good independent estimates of geopotential 
coefficients. 


APPENDIX 


ERROR ANALYSIS 


The ORAN (Orbital Analysis) computer program used for this study was 
designed for computing the effects of random and systematic errors on minimum 
variance orbit determinations. Systematic errors can be in the form of either 
adjusted or unadjusted parameters, with the effects of the latter broken down into 
effects of the individual error sources. The program computes the effects of the 
unadjusted parameters on both the recovered parameters and the orbit, with the 
orbital effects propagated from epoch to any desired prediction time. 

The program is configured for multiple arcs, with some error model param- 
eters such as station positions constrained to be common to all arcs, and other 
parameters, such as measurement biases, which differ from arc to arc. 

Force model errors can arise from uncertainties in geopotential coefficients 
through degree and order 20. Uncertainties in up to 44 individual coefficients 
can be carried, and any of these may be either adjusted or their unadjusted effects 
propagated. Alternately, or in addition, the force model error can be carried as 
the differences between complete gravity models in which case the restriction to 
44 parameters does not apply. The SAO, APL, and NWL models are built into the 
program and the differences between any two of these three, or any complete 
model supplied as input, are available as force model errors. Note that the 
gravity model difference is treated as a single parameter, and 43 geopotential 
parameters may also be considered as adjustable. Of course, adjusting a geo- 
potential coefficient removes it from the model difference set. 

Mathematically, the unmodeled error propagation is based on the following 
observations. The minimum variance orbit determination uses the basic equation 


SO - AS a + e 


( 1 ) 


to relate discrepancies (SO ) between measured and calculated observations to 
discrepancies (8a) between true and a priori estimates of the set of parameters 
to be recovered. The set 8 a includes the six orbital elements but may also 
include other parameters. The matrix A is the set of partial derivatives of the 
measurements with respect to the adjustable parameters,, and e is a vector of 
measurement ’'noise." When the least squares criterion is used to solve (1) for 
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the best estimate of a, the result is 


l 6u .r (A T WA)- 1 A T WbO, 


( 2 ) 


where W is the matrix of measurement weights. For the solution to be minimum 
variance, the weight matrix must be chosen such that 


W“ l ~ E(ec l ). 


(3) 


That is, W must be the inverse of the variance covaraince matrix of measure- 
ment noise. In the normal data reduction programs, W is generally so chosen 
because it actually is measurement random error, in which case W is rather 
accurately expressed as a diagonal matrix. 

For various reasons, the set of parameters adjusted in data reduction pro- 
grams is only a subset of those parameters having some error. For example, 
our knowledge of geopotential coefficients is by no means complete. Yet a 
truncated model is always (of necessity) used, and the error in all coefficients 
used is ignored in all variance computations. Because the net effect is that e is 
not rando " yet contains definite systematic components, we can obtain a more 
accurate ^esentation of the measurement discrepancy vector by expressing 
e as 


e = Ky + 6 , (4) 

where 7 is a set of errors in parameters previously ignored, K is the matrix of 
partial derivatives of the measurements with respect to these parameters, and 
e is the vector of measurement random noise upon which W is still based. Sub- 
stitution of (4) into (1) gives 

80 = ASa + Ky + e. (5) 


If the weight matrix for the measurements is based on e and is the same as that 
used in the data reduction program, it follows that the solution for 8 a actually 
being obtained is not that given by ( 2 ), but actually is a "biased" solution given by 
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S3 S (A t WA) A T W( 80 -Ky). 


( 6 ) 


From this relation, we may obtain by differentiation the effects of ’’unit" values 
of the set of y parameters, 



= - (A T WA )“ l A t WK. 


( 7 ) 


It follows that if the matrix K can be obtained, the effects of unit values of the y 
parameters are obtained by substituting K for the SO vector used in the data 
reduction program. A priori estimates of errors in the y parameters lead to an 
estimate of the magnitudes of the effects on recovered parameters, and the 
trajectory, of each y parameter. 

Uncertainties in the y*s are generally uncorrelated. If their correlations 
are known or can otherwise be accounted for, an estimate of the total or overall 
accuracy of the orbital solution is readily obtainable. For this study, errors in 
station locations, GM , and the geopotential were considered. (See Table 3-1.) 


REFERENCES 

1. "Resonant Satellite Geodesy Study," by B. C. Douglas md M. T. Plamiter. 
TRW Systems Final Report on Contract # NAS 5-10469, December 1967. 

2. "Orbit Resonance in the Circulation Regime," by M. T. Palmiter and G. S. 
Gedeon, TRW Systems Group Tec^ 'ncal Report# 3413.5-74, April 1968. 

3. "Theory of Satellite Geodesy," by W. M. Kaula, Blaisdell, Waltham, Massa- 
chusetts, 1966. 

4. "A Perturbation Analysis of Existing Resonant Satellites," by B. C. Douglas 
and C. A. Wagner, NASA X-Document No. X-643-68-338, August 1968. 


33 



